Functional aspects of the Hardy inequality. 
Appearance of a hidden energy 
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Abstract 



, Starting with a functional difficulty appeared in the paper |14] by Vazquez and Zuazua, 

' we obtain new insights into the Hardy Inequality and the evolution problem associated to it 

by means of a reformulation of the problem. Surprisingly, the connection of the energy of the 
new formulation with the standard Hardy functional is nontrivial, due to the presence of a 
, Hardy singularity energy. This corresponds to a loss for the total energy. The problem arises 

when the equation is posed in a bounded domain, and also when posed in the whole space. 

We also consider an equivalent problem with inverse square potential on an exterior do- 
main. The extra energy term is then present as an effect that comes from infinity, a kind of 
hidden energy. In this case, in an unexpected way, this term is additive to the total energy, 
and it may even constitute the main part of it. 



1 Introduction 



In this paper we contribute some new results on the Hardy Inequality posed in a bounded domain, 
in an exterior domain of M.^ , or the whole and on the corresponding parabolic evolution. 
The motivation came from a functional difficulty we found in the work |14j . where the following 
^ [ singular evolution problem was studied: 



(1.1) 



u 

Ut = Am + c* -j— TT, X G il, t > 0, 

u{x,0) = uq{x), for a; € ri, 

u{x,t) = in ail, t > 0. 



with critical coefficient c* = {N — 2)^/4. The space dimension is > 3 and is a bounded 
domain in containing 0, or = M^. More precisely, the authors in |l4j studied the well- 
posedness and described the asymptotic behavior of (|l.ip . Moreover, they obtained improved 
Hardy inequalities and completed the study of the spectrum of the associated eigenvalue problem. 
This problem is closely connected with the Hardy inequality: 



fl.2) / \Vu\'dx>{ — ] I ^dx, 
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which is well known to hold for any (j) € C^(0). For Hardy type inequalities and related topics we 
refer to O |9l HOl [HI Il2j . Due to this connection, c*, which is the best constant in the inequality, 
is also critical for the basic theory of the evolution equation. Indeed, the usual variational theory 
applies to the subcritical cases: ut = An + cu/\x\^ with c < c^, using the standard space 
Hq{^), and a global in time solution is then produced. On the other hand, there are no positive 
solutions of the equation for c > c* (instantaneous blow-up) , [H [5] . In the critical case we still 
get existence but the functional framework changes; this case serves as an example of interesting 
functional analysis and more complex evolution. 

In order to analyze the behavior of the solutions of Problem (jl.ip in [T3], the Hardy functional 
(1.3) Inm := ^ |V0|2 dx - (^^) ^ ^ |f "^"^ ' 

is considered as the Dirichlet form naturally associated to the equation. This form is positive 
and different lower bounds have been obtained see [HE]. Note that the expression is finite for 
u G i?o(0), but it can also be finite as an improper integral for other functions having a strong 
singularity at x = 0, due to cancelations between the two terms. To take this possibility into 
account, the Hilbert space H was introduced in [H] as the completion of the C^{Q) functions 
under the norm 

(1-4) m\Hin)=Inm, ^^C^m- 

According to Section 5 of [H], this space allows us to define in a natural way a self-adjoint 
extension of the differential operator L{u) := —Au — u/\x\'^ (the Friedrichs extension) and then 
to use standard theory to generate a semigroup and describe the solutions using the spectral 
analysis. The study of the spectrum leads to an associated elliptic eigenvalue problem, the 
solution of which turns out to be a classical problem in separation of variables. 

Problem with the singularities. The separation of variables analysis produces some singular 
solutions. In particular, the maximal singularity (corresponding to the first mode of separation 
of variables) behaves like \x\~^^~'^^^'^ near x = 0, and this function is not in Hq{Q). Now, this 
solution must belong to the space H associated to the quadratic form, hence the conclusion 
H / i?o(il). We recall that this is a peculiar phenomenon of the equation with critical exponent 
c* = [N — 2)^/4. For values of c < c* the maximal singularity is still in Hq{^). 

However, there must be a gap in the argument of [14j. We have realized that with the proposed 
definition of H, there exists a problem with the solutions of the evolution problem having the 
maximal singularity. The verification is quite simple in the case where = Bi, the unit ball in 
centered at the origin. Then, the minimization problem 

llulP 

(1.5) min !f 

u£H \\u\\j^2 

has as a solution the function 

(1.6) ei(r) = r-(^-2)/2jo(zo,ir), r = 

Jo is the Bessel function with Jo(0) = 1, up to normalization and 2:9,1 denotes the first zero of Jq. 
This function plays a big role in the asymptotic behavior of general solutions of Problem (jl.ip . 
The minimum value of p.5p is 

m = 4i. 
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Moreover, the quantity /^^(ei) is well defined as a principal value. Assuming that 
(1-7) ||ei||f, = /B,(ei), 

from the definition of H, for any e > 0, we should find a Cg°-function 0, such that 

(1.8) \\ei-<P\\l<e. 
However, setting 

ei(r) — (/)(r) = \r\^^^^'^^^'^w{r), r = \x\ , 

we have 

w{r) = Jo{zo,i r) - r^^~'^^''^(f){r) , 

which is regular at r = 0, and after some straightforward calculations we find that, using as norm 
the square root of Ibx , 

(1.9) \\e,-cp\\l= [ \x\-(^-^^Vw\'dx + l [ V|x|-(^-2)-V«;2dx. 

The last integral is not zero, due to the presence of a "boundary term", when we integrate by 
parts. More precisely, 

,1,0, I f v|.|-(--) . v.^d.= JJ(0) = 

where a; at denotes the Lebesgue measure of the unit ball in M^. Thus, 

(1-11) \Vx-^\?H>^^ -^N. 



which contradicts (II. Sh . Thus, we see that ei fails to be in iJ, since it cannot be approximated 
by Co^-functions and this will happen for every function with the maximal singularity. 

Therefore, under the assumption (jl.7p . the space H seems not to be correctly defined in [H] to 
apply the rest of the theory, since there exists a problem in dealing with very singular behavior 
near x = that is not covered by approximation with infinitely smooth functions. Actually, this 
was our first impression. 



New results 

1. The examination of the difficulty shows that the proposed norm is too detailed near the 
singularity and produces a topology that is too fine to allow the convergence of <\)n € C^(il) to 

e\. By means of a transformation already proposed in [1], we obtain a more suitable norm 

1 II 

that is equivalent to on Co(f^), but is gross enough near the singular point. In this way, we 
are able to define a possibly larger closure, that we call 'H, that contains all the functions needed 
for constructing the evolution. This is done in Subsections 12.11 and 12.21 

We proceed next to re-examine the above mentioned difficulty. We will show that the spaces H 
and H are indeed the same. What is different is the norm that was implicitly assumed to be acting 

1/2 

in H for solutions that do not necessarily vanish at x = 0, which in principle seemed to be 
taken in the sense of principal value. When both terms of /q become infinite the correct definition 
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of the norm is a particular limit that we call the cutoff limit. This is explained in Subsections 
- 12. 51 where we examine the connection of the new norm with the Hardy functional; the difference 
is characterized in terms of a certain value the Hardy singularity energy (HS energy for short) 
that we precisely define. We think that the existence of the two different norms that coincide on 
Cc(ri \ {0}) is quite interesting and was unexpected for us. 

The spectral analysis of [B^ becomes rigorous in this setting, and a contraction semigroup is 
associated to the evolution problem by standard variational methods in weighted Sobolev spaces, 
up to unitary equivalence. This analysis is carefully explained in Section [3l 

2. In Section U] we discuss a result which has its own interest; the Critical Caffarelli-Kohn- 
Nirenberg Inequalities, in a bounded domain. It was well known (see [7]) that these inequalities 
are related to the Hardy inequality with c < c^,. The critical case is as expected related with 
the Hardy inequality with c = c^,. The proper functional setting that we had to consider for 
Problem (jl.ip leads us naturally to these critical inequalities. We also give the connection of this 
new space with the Sobolev space D^'2(M^), with the use of a proper transformation. Based on 
this transformation, we describe an easy to apply argument concerning the existence of uou-Hq 
minimizers, establishing that their behavior at the origin, is precisely \x\~^^~'^^^'^ . 

3. In Section [5] we explore the existence of an analogue of the Hardy singularity energy for prob- 
lems posed in exterior domains. The Kelvin transform suggests that the most natural problem 
to study is the following: 



(1.12) 



\y\-^wt{y,t) = Au;(y,t) + c*4^, y £ B^, t > 0, 

\y\ 

w{y,0) = wo{y), for y € -B|, 

w{y,t) =0 for |y| = 5, t > 0. 



where c* = [N — 2)^/4 is the critical coefficient, = M.^\Bs{0) is the standard exterior domain 
and 5 > 0. Without loss of generality we take 5 = 1. 

Problem (jl.l2p has the striking property that the Hardy functional posed in the exterior domain 
is not necessarily a positive quantity; we will show that for functions which vanish at dBg{0) and 
behaving at infinity like |?/|~^^~^^^^ it may be negative. We avoid the difficulty by basing our 
existence theory on the unitary equivalence via the Kelvin transform. Hardy inequalities on 
exterior domains were studied in the work [2]. 

The novel feature of the bounded domain, namely, the HS energy term, does exist also in the 
case of the exterior domain, but it appears at infinity. Besides, there is a big difference with 
the bounded domain case since in this case the new energy is not only additive to the total 
energy involved in the evolution; it may even represent the main part of it. The energy term 
at infinity looks to us like a "hidden" or "dark" energy. We are not in a position to make a 
physical interpretation, but we have proved the existence of a such energy and give its precise 
formula. This seems to be the first study of an evolution problem with such curious properties 
in an exterior domain. 

Equations (jl.ip and ()1.12p are linear parabolic equations. In this sense, it is expected that 
their dynamics should be trivial. However, the presence of a singular potential with an inverse 
square power singularity changes things; due to that, the first equation has interesting behavior 
at zero and the second at infinity. Moreover, we treat the limiting case of the best constant c^,, 
in the sense of the Hardy inequality. The singularity and the best constant make the dynamic of 
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these equations far from being trivial. The appearance of correcting terms in the form of hidden 
energies has been unexpected to us. We think that the question deserves further investigation. 

4. In Section [6] we consider the Hardy functional and the corresponding evolution problem posed 
on M^. The problem was studied in ^1^, see also [13]. It is known that the Hardy inequality on 
is sharp, i. e., it cannot be improved, at least up to some L^-norm. This may be also seen as 
consequence of the transformation u{x) = v{x); we can find functions, not belonging 

to L^(M^), such that the Hardy inequality holds. To overcome this difficulty, the authors in [T3] 
made use of the similarity variables. In this case we do not find any HS energy at infinity, as in 
the previous exterior domain model, see [15] . 

The main result of the section is to prove that the Hardy functional on M^, for a certain class of 
functions, maybe improved by the L^(R^)-norm, i.e. a Hardy-Poincare inequality on M^, holds 
in this sense. The idea is to use a new, proper, transformation. 

Lemma 1.1 (Hardy-Poincare inequality on M^) Let v be any function in C^{R^), N>3 
and let 

(1.13) u= \x\~^ Jo{\x\)v, 

where Jq is the Bessel function with Jo(0) = 1, up to normalization. Then, 

(1.14) / \Vu\^dx-( ^ [ ^dx> [ u^dx. 
The best constant in (I1.14p is 1 and there exists no minimizer. 

As a consequence of the above Lemma we derive a result, see Remark 16. H which has its own 
interest. The connection of transformation u{x) = |a;[^(^^^^/^ v{x) with the reduction of dimen- 
sion from to 2 is one of the leitmotifs of paper [4i| and will be used frequently in the sequel. 
This result concerns the improvement, for a certain class of functions, of the "norm" of D^'^(M^). 
Note that D"^'^(M^) is not a well defined space. 

Lemma ll. II is. as far as we know, the first result concerning improvement, with a norm, of the 
Hardy inequality on M^. The only result, up to our knowledge, that gave an improvement of 
this inequality is in the work where the authors obtained a non-standard improvement. 

Returning now to the Hardy functional, we define a new weighted space J through (I1.13j) . We 
state some properties of this space, focused on the behavior at infinity; the functions belonging in 
this space belong also to L'^(M^), so they vanish at infinity. It turns out that ^ is a well defined 
Hilbert space. The connection of this space with the Hardy functional on is nontrivial, as it 
contains the Hardy's singularity energy at the origin, but it also contains singularity energy at 
each zero Zm of the Bessel function. Transformation (I1.13P removes the singularity at infinity. 
However, the cost we have to pay, is the existence of an infinite series of singularities at Z'fji. 
Contrary to the case of the origin, the singularity energy at these points, increases the energy. 

Some further results on Hardy type inequalities and the related hidden energies are contained 
in the forthcoming work jl3j . 

2 Proper functional setting. Bounded domain case 

We start the detailed study by analyzing the case of a bounded domain Q. C M^, N >?>. 
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2.1 Transformation and definition of spaces 

The way we follow to address the difficulty mentioned in the introduction and to properly pose 
Problem p.ip is to introduce a more convenient variable by means of the formula 

(2.1) u{x) = \x\-^^-^'>/^v{x). 

We will write the transformation as n = T{v). Clearly, this is an isometry from the space 
X = L^(r2) into the space X = L'^{dfi,Q,), dfi = \x\'^~^dx. Many arguments of |.14] were also 
based on transformation ()2.ip . which was first used in [1] and then in many papers concerning 
results on Hardy's inequalities. The great advantage of this formula is that it simplifies 
at least for smooth functions, into 

(2.2) h{v) := [ \x\-'^^-^^\Vv\^dx. 

Jn 

It is easy checked that In{u) = Ii{v) for functions u G Cq°{Q). However, the equivalence fails 
for functions with a singularity like \x\~^^~'^^^'^ at the origin, as we have hinted before and will 
explain below in detail. Our proposal is to use this formulation for the definition of the new 
space, 7i. An important observation is that when u{x,t) is a solution of equation (II. ID . then v 
satisfies the following associated equation 

(2.3) vt = \x\''-^V-[\x\-^''-^^Vv), 

with clear equivalence for x ^ 0. This last form gives the clue to the proper variational formulation 
to be followed here. First, the space associated to this equation through the quadratic form (]2.2p . 
is defined as the weighted space Ti = WQ''^{dfi, Q,), which is the completion of the Cq°{Q,) functions 
under the norm 

(2.4) ||^;||2,= / |x|-(^-2)|V7;pdx. 

^ Jn 

Following the usual procedure of the Calculus of Variations, we take an appropriate base space 
which is X = Lp'{d^,VL), and then the quadratic form ()2.2p has as form domain the subspace 
j-L where Ii{v) is finite. Then, it can be proved that L{v) := — V • (|x|~(^~^)Vf) is a 
positive self-adjoint operator in the space D{L) = {v (z Ti : L[v) G X}. It is also known that 
D^L^/"^) = j-L. Hence, the variational approach works for v. See further analysis below. 

We translate these results to the original framework. Ti is defined as the isometric space of 
Ti = M^Q'^(|x|^^^^^^dx, il) under the transformation T given by (12. ip . In other words, Ti is 
defined as the completion of the set 

^u=\x\-^v, V ec^in)^ = T{c^{n)), 

under the norm N{u) = defined by 

(2.5) l|n||2,= / [xr(^-2)|V(|x|'^tx)j2dx. 
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2.2 On the new space 

From the definition it follows that % D -?/g(r2). To show that i?Q(r2) ^ T-L, we consider any 
function in % such that u{x) ~ |x| ~ at the origin; it cannot belong to Hq{Q). We continue 
with two rather simple observations before examining the important question of relation of the 
present approach with the old one, where we will obtain a striking correction term. 

• The functions that behave at the origin hke |x|-(^-2)/2 are not the most singular in Ti. This 
was mentioned in passing in |14j . We give a counterexample here. Consider a function w which 
behave near the origin like 

v{x) {logl/\x\f, 0<a<l/2. 

It is easy to check that v G H. Therefore, belongs to 'H(ri), even if v{x) oo as 

x-^0. 

• There exist functions that belong in ng<2 but they do not belong to %. For example, 
a function v that behaves near the origin like 

v{x) ~ (log(l/r))i/2, 

fails to be in 'H, hence 

u~r 2 (log(l/r))^''^, 
does not belong to %. On the other hand, u G W^''^{^), for all 1 < g < 2. 



2.3 Connection of space "H with the Hardy functional 

By Hardy functional we refer to Iq.{u) defined in (II. 3p with the integral defined in the sense of 
principal value around the origin when both separate integrals diverge. Denote by the ball 
centered at the origin with radius e and by B'^ its complement in Q.. Assume now that u £ H, so 
that V = |xK^~2)/2 u € j-L. Then, we have that 



(2.6) IbM= I \^u?dx-(-^] I y^dx. 



By change of variables and integration by parts the following remarkable formula is obtained: 

(2.7) IbM= I \Vv\^\xf-''dx + ^^e~^''-'^ [ v^dS, 
where dS is the surface measure. Next, we denote by the quantity: 

A,(n) = ^e-(^-i) / v'dS 

(2.8) = ^e-' [ u'dS, 

that represents a kind of Hardy energy at the singularity. Is clear that 

lim / iVvP \x\'^'^dx - "-"2 



e^O « 
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In order to take the limit e — t- 0, in ()2.6p we distinguish the following cases: 

• If n € Hq{Q), then u G T-[ and we have 

A(u) := lim Ae(n) = 0, 

e-s>0 

thus the limit as e ^ 0, in ()2.6p . gives the well known formula 

In[u] = \\v\\l :=N\u), 

which holds for any u € H^{^). Note that the converse is not true; If A(ti) = 0, it does not imply 
that u € i^Q (il). For example, take a function u such that v behaves at zero like (— log 

• 11 V € T-L is such that lim|^.|_^o '^'^(^) = 'u^(O) exists as a real positive number; then u & H but 
u ^ Hq{Q,). In this case 

. X A / N N(N-2) o, , 

(2.9) K{u) = '- CON v\0), 

where ujn denotes the Lebesgue measure of the unit ball in M^. A(n) is then a well defined 
positive number and (j2.6p implies that 

(2.10) /n[u] = ||r;||| + A(n). 

We note that this is the case of the principal eigenfunction of ()3.4p and the case of the minimizer 
of the Improved Hardy-Sobolev inequality, see [15], in the radial case. As it will be clear in 
Subsection 14.11 this is the case for the minimizers of 



\\u\\h 2iV 
(2.11) min , l<p< 



u£H \ \u\\lp N - 2 

Note also that in this case the inner product of H is given by 



f ^ ^ , fN -2Y f uiU2 , 
<ui,U2>n = / \/ui ■\/u2dx - I — - — / „ dx 



(2.12) - ^^^^ ^K nVi{0)v2{0). 



• If V E is such that v at zero is bounded but the lima.^o v'^{x) does not exist, i. e., v oscillates 
near zero. For example, let 

f ~ sin ((— log l^l)'^) , |x[ ^ 0. 
Then, v belongs in Ti for some < a < 1/2, sou = |x|"(^"2)^ ^ tiiis 

case, the limit L{u) 

does not exist, since it oscillates, and from ()2.6p we have that the same happens to the Hardy 
functional, in the sense that 

(2.13) lhn{lBi[u] - Ae{u)) = \\v\f~. 

• li V £ Ti is such that liuix-^o v"^ [x) = oo. For example, let 

w ~ (- log |x|)", |x| 0. 
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Then, v belongs in H for some < a < 1/2, so u = \x\^^^^'^^v € H. Is clear that, A{u) = oo, 
and from ()2.6p we have that the same happens to the Hardy functional, in the sense that (|2.13p 
holds. 

Note that in all cases is a nonnegative quantity, for every e > and so is Ib^[u]. As a 
consequence, we obtain a generalized form of the Hardy inequality valid in the limiting case of 
()2.13p . when the Hardy functional is not defined or is infinite. We do not know if there is any 
physical meaning for the singularity energy we have found. It looks like an energy defect at the 
singularity. 



2.4 The spaces 71 and H are the same 

1/2 

We recall that H was introduced as the completion of the C^{Q) functions under the norm Iq . 
The proof oil-L = H relies on showing that the set Cq°(O\{0}) is a dense set in both spaces, and 
on observing that the two norms coincide on that subset. 

The following lemma follows from [71 Lemma 2.1], which studies the subcritical case. 
Lemma 2.1 Holds that C^{Vl\{Qi}) is a dense set in H . 
Proof By the definition of i?, it suffices to prove that 

c^{n) c Co°°(j^\{o})"-"'' 

Let p{t) be a cutoff function that is 1 for t > 2 and OforO<t<l,0< p{t) < 1. For a fixed 
u G C^{n) we define 

(2.14) us{x) = p{\x\/e) u{x) G Co~(f)\{0}). 



Then we have that 



(^j {Vue - Vuf dx 



1/2 / ^1/2 

(V {p{\x\/e)) u+{p- l)Vuf dx 



< (^J {Vp{\x\/e)f u^dx^ {p-lf\\/u\^di!^ , 

where = {x S il,e < |x| < 2e}. The first integral is estimated as 

f {Vp{\x\/e)f u^dx<c\\u\\lo.e-'^ f {p'fdx<ce^~^ f 



r2e 

{p'f dr 0, 



as e i 0. Also, the second integral tends to zero, since 

/ {p- lf\Vu\^dx <c dx, 

so 

(2.15) / {Vue - Vuf dx 0, 

as e I 0. In addition. 



/ 



n \x? Jb. 
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as e 4, 0. Finally, we conclude that 

\\Ue - u\\h 0, 

and the proof is completed. ■ 

Next we will prove that the CQ°(r2\{0})-functions are also dense in 7i and hence in 7i. However, 
the functions defined by ()2.14p are not good approximations inH. In other words, the arguments 
of [71 Lemma 2.1] cannot hold in the critical case a = (N — 2)/2. 

Lemma 2.2 Let v € Cq^{Q,), v{0) ^ 0. Then the v^, given by ^^(x) = p{\x\/e)v{x), cannot 
approach v in 7i. 

Proof Let p be as defined in Lemma [2?T1 Let v G C^{0,) and V£{x) = p{\x\/e) v{x). Then 
(2.16) +/ [xr(^~2)(/,-l)2|Vt;|2dx, 

where = {x ^ Vt^e < \x\ < 2e}. Integrating by parts the second integral and noting that 
p'{2) = and p'{l) = 0, we have that 

/ \xr^''~^^{p-l)p"v'dx 
-e-i / \x\-^^-^^{p-l)p'v''dx. 

Thus, we obtain that 



\\ve-v\\l = -e-^ j \x\-^''-^^ i^^p" + p'^ {p-l)v^dx 

(2.17) + / |xr(^-2)(p-l)2[V^;pdx. 

JA^ 

Taking e << 1, we have that v'^{x) ~ ^^(0), x € Af,. Then, for t = r/e, the first integral in (j2.17p 
is equal to 

1 / Lxl-(^-l) {\^n"^n' 



-e / \x 

'a. 



(^-^p" + p'^ {p-l)v^dx = c [tp'{t))' {p{t)-l) dt + 0{e). 



Integrating by parts the last integral and using again that p'(2) = p'(l) = 0, we obtain that it is 
equal to 

(2.18) -^'^ f i^^p'' + P'^ {p-^)v'^dx = c j\[p'{t)f dt + 0{e). 
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On the other, for e J, 0, the last integral of ()2.17p 

j-2e 

\x\-^^^'^) {p-lf\\Jv\^ dx <c I rdr^O. 

'B2e Jo 

Hence, from (|2.17|) and (|'2.18p we conclude that 

(2.19) Y\ui\\ve-v\\% = c tip\t)fdt, 

which is a positive quantity. The proof is thus completed. ■ 

The special cutoff functions that are dense in 71 are the ones that allow to prove that {0} has zero 
capacity in two space dimensions. The connection of transformation ()2.ip with the reduction of 
dimension from to 2 is one of the leitmotifs of paper [4J . More precisely, we prove the following: 

Lemma 2.3 The C^{Vt\{Q}) -functions are dense in T-i. 

Proof By the definition of Ti, it suffices to prove that we can approximate smooth functions 
vanishing at dil. by smooth functions that vanish also near x = 0: 



Let Pe{t), e G (0, 1), be the cutoff functions defined as follows: (i) Pe{t) = 1 for t > e; (ii) Pe{t) = 
for < t < e^; (iii) in the remaining region 1/e^ < t < 1/e, it has the special form 

Pe{t) = Ce\og{t/e^), Ce = (log(l/e))-l . 

For a fixed v € C^($7) we define fe(x) = pe{\x\) v{x). Note that 

iVpedxDI = for X e As = {x £ ^ : e'^ < \x\ < e} 

\x\ 

being zero otherwise. Then we have that 



2 



< 2 [ |xr(^-2) |Vp,(|x|)pt;2dx 



I A, 

(2.20) +2 [ \x\-'^^-^^ {1- p, f\Vv\^dx, 



Be 



We will prove that letting e J, 0, the integrals in (|2.20p go to zero. The only one that is delicate 
is the first. We have 

(2.21) [ \x\'(^-'^\Vps{\x\)\'vUx<CMl f 4- = C\\v\niog{l/e))-^ 

J A, Je^ r 

and this tends to zero as e — )• 0. Approximating the functions v^, which vanish around the origin, 
by Cg°(il\{0}) is now standard. ■ 

From Lemmas 12. II and 12. 31 we have that the spaces H and T-L may both be defined as the closure of 
C^(O\{0})-functions, with respect to (II. 4p and (12. Sp . respectively. However, for such functions 
these two norm are equal, i.e., 

(2.22) \\w\\l^ = Iu{w), 
for any w G C7~(J]\{0}). Thus, 
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Proposition 2.1 The spaces H and Ti coincide. 

For the space which is defined by ()2.4p . we have that 

Lemma 2.4 The space % contains all the functions that satisfy v\q^ = and \ < oo. 

Proof We first note that the bounded functions with vg^i = and < oo belong in T-L. This 

follows using the Vs functions defined in Lemma 12.31 see ()2.2ip . 

Next, for any v with vqq = and \ \v\\^ < oo, we will prove that v is approximated by a sequence 
of bounded functions in 7i. Indeed, if we define Vn as follows 

Vnix) = v{x), if \v{x)\ < n, Vn{x) = n, if \v{x)\ > n, 

we have that 

(2.23) / \x\-^^-^'>\V{v-Vrr)fdx= j \x\-^^-^^ \Vv\^ dx < OO. 

Now, it is clear that the sets = {x E : \v{x)\ > n} form a monotone family and the measure 
tends to zero in the limit n ^ oo. This means that the above integral goes to zero as n ^ oo 
and the proof is complete. ■ 

We have also obtained the above results following a different approach. This is included in the 
Appendix since some readers may be interested in that approach. 



2.5 On the norm of H 



Let us examine further the definition of the norm N that will be considered for the space H = T-L. 
We know that the norms I^'^ and N coincide on functions Hq{Q), and also that In{u) is larger 
than N'^{u) = when they differ. More precisely, 

(2.24) \\u\\jj = lim {Ibc[u] - A,(u)) . 

Now, if for any u (z H we consider a sequence of cutoff approximations Us{x) = pe{x)u{x) with 
as in Lemma (12. 3p . then = lQ.{ueY^'^ and tij — y u in H. The limit value 11^11^ — 

lime_^o ll'Uelll/ is what we call the cutoff value of the Hardy functional, and produces the correct 
norm on H. 

As a conclusion, the space H as it is defined in Vazquez-Zuazua [14j is a well defined space, as 
the completion of Co°(r2)-functions with respect to the norm \\4>\\\j = £ Cq^{^). In 

this space there exist "bad functions" , such that {u) defined as an improper integral does not 
coincide with the limit of the sequence of cutoff approximations. Even more, it can happen that 
the principal value of the integrals in is not well defined, either oscillating or infinite. For 
example, let u behaving at the origin as Ixl^*-^^^-*/^; for this function, the quantity Iq,{u) is well 
defined, but its norm in H is not but it is equal to Iq.{u) — A(n). For the first eigenfunction 

1/2 

eii the norm is not (ei), but 

lleilll^ = /^(ei) - A(ei) = /n(ei) - ^^^^ — '^^n = fJ-i- 
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As a result the minimization problem (jl.Sp and the following one 

(2.25 mm-- — -rrr- 

u&H ||m||22 

are not the same, since the first admits a minimizer while the other does not. This will also happen 
in the case of the minimizers of (|2.1ip , see Subsection 14. 1|, and in the case of the minimizer of 
the Improved Hardy-Sobolev inequality in the radial case, |15j . 

As it is stated in |14l pg. 127], the solutions of the stationary and evolution problem in the 
critical case, c* = {N — 2)^/4, can be obtained as the limit when c t c* of the solutions of 
the subcritical case c < c*. This is true in the space with norm (12.24p . and not through 
I. For example, let ei^c be the eigenfunctions corresponding to the principal eigenvalues of the 
subcritical case, on the unit ball. In |l4i Theorem 3.3] the exact formula of these eigenfunctions 
were obtained: 

ei,c(r)=r-(^-2)/2j^(z^,ir), 
where m? = c* — c > and Zm,i is the first zero of the Bessel function Jm- Then 

llei.clln = In{ei,c{r)) /^(ei) - A(ei) = ||ei||^, 

as c '\ . 

This strange situation has to do with the definition of In{u); the existence of the difference of 
two improper (and comparable) integrals. This is the difference between the case of a classical 
Sobolev space or the case of the Hardy functional with constant less than the best constant. 



3 Application to the evolution problem 

We now justify that the results described in [ll] for the solutions of Problem (jl.ip hold in the 
space Ti defined in (j2.5p . with Q a bounded domain of M^, > 3. Is clear that this space is 
a Hilbert space and, as stated in [HI Theorem 2.2], T-L is imbedded continuously in the Sobolev 
space VF^'<?(17), 1 < g < 2, i.e., 

(3.1) \\u\\n>C{q,n)\\u\\wi,q(^Q), l<q<2. 
Thus, the compact imbedding 

(3.2) n^LPin), i<p< 



A - 2 

holds. Moreover, we may justify that all the results concerning the spectrum of the related 
eigenvalue problem given in [14j hold for H. For this it is sufficient to proceed as follows. Assume 
the Cauchy-Dirichlet problem for equation 



(3.3) ut = Au+ 

and the corresponding eigenvalue problem 



A-2\^ u 



(3.4) Au+{^—]--^ + nu = 



/ N-2 y u 
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Setting It = l^l 2 V, (|3.3p becomes equation (|2.3p away from zero. The distributional operator 
L given by 

(3.5) L(t;) = -|x|^-2V-(|xr(^-2)V7;) 

is an isomorphism from % into its dual and from (j3.2p we have the triplet 

Then, by restriction, we may define the surjective operator : D{L^) d X ^ X , with domain 

= {/ G ?^ : 1x^-2 V • (|xr(^-2) V/) G x] . 

If we denote with L, is easy to see that L is self-adjoint with compact inverse. So, we can 
form an orthonormal basis of eigenfunctions in X, with eigenvalue sequence 

(3.6) < /ii < /Z2 < ... < /i/fc < ... oo. 

(As in [H], separation of variables allows to find explicit formula for the eigenfunctions and 
eigenvalues when i7 is a ball). Returning now to Problem (j3.3p . we have that the result of [14| . 
concerning the spectrum of the associated eigenvalue problem as well as the asymptotic behavior 
of the evolution problem, both are true. 



3.1 Energy calculations 

In order to justify the above choice we make the following calculations. Let be the eigenfunc- 
tions corresponding to (j3.6p . Assume that v{t^x) is a solution of (I23D and v G L2(0,T : U). If 
f(t,x) = Cfc(t) we have that 



Hm\ = Y.^lit) and \\v{t)\\\ = Y,cl{t)\m\\~. 



I ~ 

k=l k=l 



If we denote by E[t) := \\v{t)\\^, we obtain that 



(3.7) ^ = 2f;c,(t)4(t). 

k=l 

and for E^^^{t) :=< v{t),(j)k >jj we get that 



(3.8) ^^ = 4(t). 
However, the evolution equation implies that 

(3.9) ^^^ = -<vAk>^=-ckm\i. 



Thus, dSSI) and (I3J|) imply that c'^{t) = -Ck\\(t)k\\\- Then, from ([321) we conclude that 
(3.10) ^ = -2 f; cl{t) = -2\\v\\l = -2 \Vv? dx. 



k=l 
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• Let us now see where the original Hardy functional fails by performing this energy calculation 
in terms of u. To avoid the singularity we calculate the energy in il.^ = il\i?£(0)with e > small. 

d f 2 , /■ A , /■ , 

— / u ax = I u /\u ax + / - — ax 

,2 



/ (iVnp — c*-j — ) dx + uUndS. 
Jn, V kr/ J\x\=e 



'\x\=£ 

Using u = V r^^^'^'^^^'^n with r = \x\ > 0, we calculate the last integral 

uun dS = -!^ei-" I v'^dS + e^-" f WrdS. 

r=£ J r=l 



2 

Recalling now formula ()2.7p we arrive at 

d 



^,j^-'dx = -\\v\\l^^^^^+e-~- j^^vv..dS 

If we now pass to the limit e — >■ and realize that the last integral goes to zero whenever v (zTi, 
we get formally 

This is precisely the choice made by our setting. Let us briefly justify the term that went to zero 
in average when v is bounded: integrating in e from S to 26 and putting As = {x : 5 < \x\ < 26}, 
we have the following estimate for the average 

1 /■2<5 I- (-If- r / r \^/'^ / r'^^ \ ^^"^ 

\-J e^-'^ds J vVrdS\<jJ \Vv\\x\^-''dx<j i^J \Vv\^ {xl^-^'dxj i^J rdrj 

that goes to zero as 5 — >• 0. So, the weak formulation (or the energy equation) of (jl.ip is the 
following 

^juf = -||n||^ = ~ J™i i^Bsiu] - ^eiu)) 
for every u & H. The space H is really the energetic space. 



4 Further properties of the spaces 
4.1 Further properties of T-L and "H 

We observe that imbedding (|3.ip gives the following corollary, which completes the results ob- 
tained in [6] (see also [7J) concerning the Caffarelli-Kohn-Nirenberg Inequalities, in a bounded 
domain, in the limiting case where a = ^^y^- 

Corollary 4.1 (Critical Caffarelli-Kohn-Nirenberg Inequalities) Assume that is a bounded 
sequence in Ti. Then Un = \x\~^^~'^^^'^ Vn is a bounded sequence in %. The compact imbeddings 
(|3.2p imply that, up to some subsequence, converges in L'^iVt) to some u. Thus, we obtain the 
compact imbeddings 

~ 2 /V 

(4.1) LWx\-'i'^^-^'^/'^dx,^l), for any I < q < ——. 
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Then, for every < s < ^^^-2^q, we further obtain the compact imbeddings 

2N 

(4.2) L'ii{\x\-'dx,n), for any l<q< 



N -2' 



where the weighted space L'*{w{x) dx^Vl) is defined as the closure ofC^{^l) functions under the 
norm 



Li{w{x),n) 



1 

w(x) \(t>\'^ dx 

n 



Remark 4.1 In (14. Ih it is clear that q cannot reach jf^- For this value of q, the best that we 
can have are Improved Hardy- Sobolev Inequalities, see and the references therein. 



defined as the closure of C^(M^) functions under tlie norm 



In addition, we can relate these spaces, in the radial case, with the space I'^'^(]R^), which is 

roc 




^i.2(iRiV) = / |V0pdx. 

If we denote by Hri^), and Vr''^{M.^) the subspaces of %, H and 2?^'^(]R^), respectively, 
which consist of radial functions, we have that 

Proposition 4.1 For some function v € T-LriBR) we set 

-log (^y 

Then, v G HriBji) if and only if vu € Pr'^(IR^) and 
(4.4) \\v\\^^^^^^ = {N-2)-'\M^r,. 

Observe that (I4.4p is independent of the radius R and in the case where = 3 the norm in 
T-Lr{BR) coincides with the norm in Moreover, the definition of Ti and ()4.4p imply that 



Corollary 4.2 For some function u we set 



1 

x\\\ 



(4.5) w{t) = \x\ 2 u{\x\), t=j^-log(^y 

Then, u € TiriBji) if and only if w ^ Pr'^(M^) and 

(4.6) IHlL(B«) = (A^-2)-MHlii,2(K.). 

Transformation ()4.3p was used in |15j . For a discussion concerning the construction of such 
transformations for general Hardy inequalities we refer to [l3j . 
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4.2 Nonexistence of i^Q-mininiizers 

Transformation (|4.5p provides us with an extra argument concerning the nonexistence of Hq- 
minimizers. And not only this; we are able to obtain the exact behavior of these minimizers at 
the singularity. We will prove that these minimizers belong to H, they do not belong to Hq and 
their behavior at the origin is exactly 

As an example, we will prove that the minimizer of (|1.5|) cannot be an i^Q-function. As already 
mentioned in the introduction, this minimizer is the function ei, given by ()1.6p . and does not 
belong to Hq. 

Assume to the contrary that u £ Hq is a minimizer of p.5|) . Then, u may be chosen to be a 
nonnegative and radial function, i.e. u{x) = u{r) > 0. Let w be the transformation of u, through 
(14.5(1 . Since u G Hq we obtain that 
(4.7) w{0) = 0. 

Moreover, we have that w € D^''^(R^) is a minimizer of 
(4^g) 1 IuN\Vw\^dx 



(AT -2)2 f^,V{\x\)w^dx' 

where l^(|x|) = |rc|^2(^-i) e-2|a;| Note that if we set V{0) = 0, y is a continuous function. 

Then, w should be a nonnegative solution of the Euler-Lagrange equation corresponding to 



-Aw = c{N)V{\x\)w, w£D^'^(R^). 

However, the application of the Maximum Principle contradicts (|4.7p . hence (|1.5p cannot have a 
iifQ-minimizer. 

This argument may be applied to more general problems; 



Proposition 4.2 Let ^ he a hounded domain of M^, > 3, containing the origin. Then, 
minimizers of 

(4.9) mm— \<q< 



cannot exist in Hq{Q,). 

Proof The compact imbeddings (|4.2p imply that the minimization problems (j4.9p have a solution, 
let us denote it by CQ^q, and this solution is attained by some function 

If il. is radially symmetric, then the conclusion follows using the same argument that was applied 
for ei. In the case where is an arbitrary domain, we make a symmetrization that replaces f2 
by a ball Br with the same volume and the function UQ^q by its symmetric rearrangement ^. 

Assume that UQ^q G Hq{Q,), then = Iniun^q]. It is well known that the rearrangement 

does not change the L'^-norm, increases the integral J v? /\x'\^ dx. Also, decreases the Hq{^)- 
norm, so that ^ belongs to Hq{Br). Finally, we get the contradiction that Uq ^ is the minimizer 
^*_Bfl,(j and belongs to Hq{Br). Thus, u^j^g cannot exist in HQ{yL). ■ 

The case q = has the same quantitative behavior (in the radial case) and this maybe 

obtained following the same argument. 
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As mentioned before, the case of ei and the minimizer of the improved Hardy-Sobolev Inequahty 
(in the radial case), belong to the second case of Subsection 12.31 This means that they behave 
at the origin like \x\~^^~'^^/'^ . Thus, the Hardy functional for these functions is a well defined 
positive number, although it does not represent their //-norm. These functions do not belong 
to the "worst" cases, where In is not well defined or is infinite. As a corollary of the previous 
argument we have that the same happens to every minimizer uq^„ of ()4.9p . 



Corollary 4.3 Every minimizer u^^q of (j4.9|) behaves at the origin like \x\ 

Proof Assume the opposite. Since does not belong in H^, in the sense of Subsection 12.31 we 
have that the Hardy functional of UQ,^q is not well defined or is infinite: 

lBA^^,q] oscillates or is infinite, as e — )• 0. 

Since, uq^^ € H{Vt), from ()2.24p we have that the same happens to Ae(nQ^g), as e ^ 0. 

We consider first the radial case, i.e. un,q is a radial function and VL = Br. Let wn^q be the 
transformation of UQ^q, through (j4.5p . Then, from the definition (j2.8p of A^, WQ,g(0) is not well 
defined or is infinite, respectively. On the other, wn^q € -D^'^(M^) is a minimizer of 

1 f^N |Vwp dx 

(^-^y^'^' {f^.Vq{\x\)m^dxf'' 

i.e., wn^q is the solution of the corresponding Euler-Lagrange equation 

-Aw = c{q, N) Vq{\x\) \w\'^~^ w, w e D'^'^{R^), 

with Vq{\x\) = |x|~2(7V-i) j(^) ^ N-{N-2)q/2 > 0. Vq is a continuous function, 

if we set V^(0) = 0. However, standard regularity theory implies that this is a contradiction. 
Next we treat the general case. We consider Uq the symmetric rearrangement of u^^q. In the 
case where lima;_j.o |x|^^~^^/^ UQ^q{x) = oo, ^ does the same at zero, since their level sets have 
the same volume. In the case, where lim^^^o \x\^^~'^^^'^ UQ^q{x) oscillates, we make the unproved 
hypothesis that Uq ^ cannot be defined at 0. However, Uq ^ should be a radial minimizer and this 
is a contradiction. ■ 

As in the case of ei , we emphasize the fact that the minimization problem (j4.9p and the following 



mm ■ 



UGH |L,||2/9 

are not the same. 



5 The case of the exterior domain 

We consider Problem ()1.12p describing the evolution (up to some weight) of the Hardy potential 
in an exterior domain. We may fix 5 = 1. Our arguments will be based on the unitary equivalence 
with the previous problem posed on a ball. For that we use the Kelvin transform in the form 

(5.1) u{x) = \y\^'^ w{y), y = . 
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These formulas transform solutions u{x,t) of Problem (jl.ip . i.e. 
(5.2) 



1l( X t] 

ut{x,t) = Au(x,t)+c^ , , X £ ^i(O), t > 0, 

\xr 



u{x,0) = uo{x), for X E Bi{0), 

u{x,t) = in dBi{0), t>0. 

defined in the unit sphere -Bi(O) into solutions w{y,t) of Problem (I1.12j) posed in B^{Q), its 
complement in M^. We will write the transformation as u = K,{w). Note that for smooth 
functions holds that 

A.u(x) = |yr+^ i^yw{y) and % = \y\^^^ 

fI \y\ 

The equivalence of the equations is clear equivalence for x 7^ 0. The differences will appear in 
the energy near the singularity versus the "energy at infinity". 

Basic properties. We will address the questions of existence of solutions for Problem (jl.l2p by 
means of this equivalence, that will also be used as a clue to the proper variational formulation 
that will be followed. As we saw, the functional space which corresponds to (|5.2p is H with norm 
given in ()2.24p . Our proposal is to use this formulation for the definition of the new space, W. 
The space W is defined as the isometric space of H under the Kelvin transformation (j5.ip . In 
other words, W is defined as the completion of the set 

w^(y) = |yr^+'n(^), uec^iBM). |y|>i 

under the norm HtfUw defined by 

(5.3) ||u;||vv = lim (/bj(o)\b,M - Ae(n)) , u = lC{w). 

The first eigenpair of the corresponding eigenvalue problem 

(5.4) — Aw — c*7-^ = lyl""^ 

\y\ 



(5.5) /XI =4^, ei = M-(^-2)/2jo 



IS 

.2. ^, = U,|-(^-2)/2 ( 

\y\ 

We have the following Result: The well posedness of p.l2p in the space W is understood through 
the unitary equivalence with H. The existence, uniqueness and stabilization results of Problem 
K2h apply for Problem (frT2]l . 

Hardy functional. Next, we investigate the connection of the space W with the Hardy func- 
tional, /_Bj(o) defined as: 

,2 . ^ 



(5.6) I (o)[</,]=/ \y^l^dx-[^— / f-^dx, 

which is positive for any compactly supported (p & C°°(i?f(0)) that vanishes on the boundary. 
We denote by Is and by Ii/e the Hardy functional defined on Bi{0)\Bs and Bi/s{0)\Bi{0), 
respectively. Recall that Ae(u) is given by ()2.8p . 
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Lemma 5.1 We have the following fundamental relation: 
(5.7) /,[n] = /i/JH + 2Ai/.M, 

where 



N -2 
^i/eiw) = — — e 



dS. 



Moreover, it is clear that if u = IC{w) then 

(5.8) Ae{u)=A,/,{w), 

Proof. In order to justify Equations (|5.7p and (j5.8p we proceed as follows. Let w € W, i.e. there 
exists u (z H, such that u = JC{w). Then, 



/ \Vxufdx 



-2N+4 



Bi(0)\Be(0) 



dx 



+{-N + 2f [ 

(5.9) +{-N + 2) [ 

Integrating by parts the last integral we get that 



-2N+3 ^ . ^ 2 



(0)\B.(0) 



X 



dx. 



\-2N+3 



(0)\BaO) 



^.V,w^-^] dx = (iV-2) 



I I-2N+2 2 
X\ W 



-2N+3 



Bi(0)\i?,(0) 

2 f X 
W —rr 



dx 



dS. 



Then from (15.91) we have that 



IVrUh dx 



Bl(0)\Be(0) 



-2N + 4: 



Bi(0)\B,(0) 



VxW 



(5.10) 



X 

2 ' ^ 



+(iV- 2)5-^^+3 ^ w'(j^) 

JdBJO) J 



dx 
dS. 



However, the determinant of the Jacobian of the Kelvin transformation in d > 2 dimension is 
equal to — jxp*^ and moreover 



X 



\X\ 



\x\-^\Vyw{y)\\ 



Finally, we conclude that 



IVrUl^ dx 



Bi{0)\B,(0) 



Ibi 

+(iV-2)e / 
JdB. 



/.{0)\Bi(0) 



\Vyw{y)\'^dy 



w\y)dS. 
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Thus, dSZD and 1^ hold. □ 

When we apply these results to functions in the class W (by density), we are able to give the 
following unexpected definition of the norm of W, 

(5.11) llu-ll^ = lhn(/i/,H +Ai/,(u;)) . 

So, the weak formulation (or the energy equation) of (jl.l2p translates into 



1 

- / Wf = —\\w 



,,2 



VV = - 1™ {h/eH + ^l/eiw)) , 



for every w (zW. 



Conclusions and remarks, (i) We have shown that a correcting term also appears as in the 
energy analysis of the problem posed in the exterior domain. Actually, the correcting term has 
the same absolute value as the Hardy singularity energy considered in the problem in a bounded 
domain, but now it represents a kind of energy at infinity. However, there is a big difference from 
the bounded domain case since in this case the singular energy acts in an additive way to the 
usual Hardy integral. 

(ii) Moreover, as we will see below, this new term may be the main part of the total energy, since 

may be also a negative quantity: we do the calculations for ei, given in (15. 5p . 

Il^lllvv = = m = ^0,1 - ^•'^6 

However, 

lleillw = 4jv\iJi(o)(ei) + A(ei). 

Thus, 

N{N - 2) 

T 

which is negative for = 3. 

(iii) As we saw in Subsection 12.31 A^/^ may have a bad behavior; oscillating or tending to infinity. 
In this cases, the Hardy functional Ii/^ becomes negative and has the same behavior with Ai/^, 
so that the sum of them to become a positive real number. 

(v) As a conclusion, we can say that A is a hidden energy that "comes" from infinity and is not 
only a gain of the total energy but it may represent the main part of the total energy. We are 
not in the position to say that, but we feel that this energy may represent a simple form of "dark 
energy" . 

(vi) Assume (jl.ip posed on whole M^. Then, the HS energy is present only at zero and not at 
infinity, cf. next section and [H] and [T3]. However, the Kelvin transformation implies that the 
corresponding problem for Equation ()1.12p is well defined on M^, and it can be checked that the 
HS energy is present as a dark energy at infinity, and not at the origin. For the details, we refer 
to [13]. 



%Jv\Bi{o)(ei) = 5.76 UJN, 
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6 A Hardy-Poincare inequality on 



Throughout this section we refer to Hardy functional I^n , defined on M , > 3, 

(6.1) Irn[u]= \Vu\^dx-l^—] / —dx, 

Jm V 2 y j^N \x\^ 

with limits taken in the sense of principal value if both integrals diverge. First, we discuss the 
reason why the transformation ()2.ip is not proper in the case of M^. Let u E C^(M^) and v 
given by (j2.ip . Then this formula simplifies I^n{u) since: 

(6.2) I^^(u) = Ii^n{v)= I \x\-^^-^'^\Vv\'^dx. 
However, for radial functions Ii^.'^iv) is equal, up to some constant, to 



and it is well known that L'^'^(M^) is a not well defined Hilbert space; it contains the constant 
functions, so we cannot have an imbedding in L'^{M?), for any 1 < q. Thus, defining the space 
?^(M'^), as in the bounded domain case through Ii^^n{v), we get that this is not well defined 
since it contains the constant functions. This also a way to see that the Hardy functional on 
cannot be improved by the L^(M^) norm, for any Cq°(M^)- function. 

Definition of the new proper transformation As it will be clear, the transformation (|2.ip 
is not the only proper one to simplify the Hardy functional. The way we follow to solve the above 
difficulty is to introduce another convenient variable by means of the formula 

N — 2 

(6.3) u{x) = \x\ 2- Jo(|j;|) v{x) := Tm,n{v), 

where Jq is the Bessel function. Note that Jo(0) = 1, up to normalization and Jq ~ r"^/^ cos(r — 
|), for large r. 

We will write the transformation as u = Tg^N{v). Clearly, this is an isometry from the space 
X = L2(M^) into the space X = L'^{dn,R^)), dfi = j2 dx. The great advantage of this 

formula is that it simplifies I^n{u), at least for smooth functions; For any v € C^(R^) holds 
that 



4^[n] = / \x\'^^-^U^{\x\)\Vv\^dx 
(6.4) + / \x\-'^^-^^ J^{\x\)v^dx + A{v), 



where A is the Hardy singularity energy be defined in (|2.9p . 
Proof of Lemma [777] Let f be a C^(IR^)- function, and set 

(6.5) u = \x\ 2 Jo(|a;|) v, 

where Jq is the Bessel function. If we do the change of variables (|6.5p . we have that Jqv €z 
C^(M^), so from Subsection O 

(6.6) I^n[u] = [ |:rr(^-2) {V {M\x\) v) f dx + A{v). 
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After some calculations, we get that 

I^n[u]= [ \x\-'^^-^^ J^{\x\)\Vv\'^dx+ [ (VJo([x|))2?;2(ix 

Jr!^ JrN 

(6.7) +[ \x\-^^-^^ Jo{\x\)VJo{\x\)-Vv^dx + A{v). 

If we denote by to be the domain M^\{B,(0) U Br{0)}, holds that 

Jq{\x\)VJo{\x\) -Vv^dx = - [ V (\x\~'^^~^^ Jo{\x\)VJo{\x\)) v^dx 

+ / \x\-^^-^^ Jq{\x\)^ Jq(\x\)v^ r]ds 



Taking the limit as i? ^ oo we have that the second integral tends to zero, since v € C(5"(R^). 
For the last integral, we have that Jq = —J\. Where J\ satisfies Ji(0) = 0. Hence in the limit 
e J, 0, the last integral also vanishes. Then, from (|6.7|) we have that 

= / \x\-^^-'^^ Jl{\x\)\Vvfdx 

Jo([x|) f AJo(|x|) - ^^J^(lxj)) v^dx + A{v). 



However, 

iV - 2 

AJo(r) Joi''') = —Joi'''), for any r. 

If we denote by J^n [v] the quantity 

f \x\-^^-^^ J^{\x\)\Vv\'^dx+ [ \x\-'^^-^Ul{\x\)v^dx, 

we conclude that the Hardy functional is equal to 

(6.8) I^n[u] = J^n[v] + A{v), 

thus 

/ |Vn|2dx-f^^V / ^dx > [ J2([x|)r;2dx 

JrN \ 2 J J^N J^N 



dx. 

Jrn 

Finally, we have that 

. ^ r Irn\x\-^^-^^ Ji{\x\)\Vv\^dx 

To see this we use the standard procedure; Let Vn € Co(M^) be such that Vn{r) = Ci for 
< r < , \v'J < Cg for ^ < r < ^^^^ and t;„(r) = for r > (Jl:^. Then, using the fact 
that Jo{r) ~ r^-*^/^ cos(r — j), for r >> 1, we have that 

^ (n + l)7r 



Ixl-^^-^') J^{\x\)\Vvn\^ dx < c I * cos^ir - ^)dr = C, 



4 
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for any n, large enough. However, 

(n+l)7i 



|x| Jq{\x\) v"^ dx > c I cos^(r — — ) dr ^ oo, 







as n ^ oo. Hence, 1 is the best constant in p.l4p and there exists no minimizer. ■ 

Remark 6.1 The above arguments may be applied to the case of D^''^(M?): Let v he any function 
in C^(M2) and 

(6.10) u = U\x\)v, 
Then, 

(6.11) / \Vu\^dx > / u^dx. 



The best constant in (|6.11|) is 1 and there exists no minimizer. 

Remark 6.2 In order to clarify the transformation (j6.3p we argue as follows. Let C be a 
hounded smooth domain. Assume that (Ai,ni) be the principal eigenpair of the Laplace operator 
on Vt. Then, for every v € Cg°(i7) and u = uiv, holds that 

||n||^i|.^-, = / |V?j| dx + Xi / Ui\v\ dx. 

In our case, the Bessel function Jq solves the following eigenvalue problem 

~Au = \u, xGM^, 

with A = 1. 



The following remark, makes clear that the transformation (j6.3p is not valid for every u € 
C(^(R^), so that J[v] to be finite. Thus the Hardy-Poincare inequality (I1.14p . in Lemma ll.l^ 
makes sense for a certain class of functions. These functions, as it will be clear in the sequel, 
behave at the origin like the functions of H, at dBz^(0), the zeroes of the Bessel function, must 
vanish like ||x| — z^l", a > 1/2, and at infinity must tend to zero. 

Remark 6.3 The change of variables id. 5\) is not proper for every u € Cg°(M^). For example, 

JV-2 _-i 

if supp{u{r)} E (0,zo,i), then v = \x\ 2 {\x\)u has J[v] < 00. While, for supp{u{r)} = 
[e, zi + 5], Q < 5 < Z2 — zi, such that u{zi) ^ 0, we have that J[v] = 00. 



Functional spaces Next, we discuss some properties of the corresponding functional space. 
We focus on the behavior at infinity. We introduce the space = VF^'^(ci/i, M^), d^ = 
|^|2-Ar jS^'i^.j^ ^j.^ ]-,g defined as the closure of C^(M^) with respect to the norm 

||^;||2 := jj^^[„]= / \xr^^-^U^{\x\)\Vv\^dx+ [ Ixl-^""-^^ J^(\x\) dx. 

We can easily verify that this weighted space JT" is reasonably defined (see, for instance p^), 
since the weight Jo(|x|) belongs to Lj^^, with Jq"^(|2;|) belonging also to Lj^^. 
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The latter is true since the zeroes of the Bessel functions are simple. Moreover, as we will see 
below, every function belonging in J tends to zero at infinity. Thus, J" is a well defined space. 

We examine the behavior of the elements of J at infinity. Note that close at the origin, J is 
equivalent to H. 

Radial Case Let A = {R, +oo)\ |J(2;m — e, + e) C M, e small enough. Assume that v = v{r) 
is a radial function in J and R is large enough. Then, 



poo 




cos^(r - 7r/4) ((v')^ + f^) dr. 



This means that in -B^\IJ (i?2„+£(0)/i?z^_£(0)) the norm of v is equivalent to ||f From 
this we conclude that 

lim v{r) = 0, 

r— >oo 

up to the behavior at r = z^- Moreover, 
for any 2 < q < oo. 

Nonradial Case Let w be a zero mean value function in J and u = T^n{v). Then, (|6.8|) holds 
with A{u) = 0. Since u is nonradial, such that /Kiv[ti] < oo, we may obtain that u G D^'^(R^), 
thus u € (R^). Moreover, from Lemma [1.1 1 we have that u belongs also to L^(M^). Finally, 
we get that u = o(|x|~^/^) and so v = o{\x\~^^'^), as |x| — >■ oo. 

General Case The previous discussion was made in order to highlight some further properties 
of i7. We can argue directly from the definition of the norm of J^. Finally, 

Lemma 6.1 Assume that v ^ J . Then, 

lim v{x) = 0, 

up to the behavior at dB^^ (0) . 

We translate these results to the original framework. J is defined as the isometric space of J 
under the transformation T^n given by (j6.3p . The above definition is equivalent to say that J is 
defined as the completion of the set 

|n = \x\-^ M\x\)v, V € Co^(M^)} = rM^(Co°°(M^)), 
under the norm N^n{u) = \\u\\j defined by 

\\u\\\ =: ( |xr(^-2) j2|V(|x|^ Jo-in)pdx 
(6.12) +/ J-^updx. 

Note that the last integral is nothing else than the L^(R^)-norm. Thus, for any u € ^7, we have 
that u decays to zero at infinity. 



The case of J We discuss the connection of J with the Hardy functional. We also give the 
properties of the functions belonging to J . 
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Lemma 6.2 Let i?£,m = -^^(O) U^^^^ i?2,^+e(0)/i?2^_£(0), Zm are the zeroes of the Bessel func- 
tion. For any u ^ J , holds that 

oo 

(6.13) ||n||^ = limlijc^^[?i] - Ke{u) + ^ K+,m{u) - J^e-,m{u), 

m=l 

where, ^ denotes the complement of B^^m in , is the Hardy singularity energy at the 
origin, defined in (|2.8|) . and 

K+,m{u)=l jQ^J^U^ds, K-,ra{u)=l Jq^ J'qU^ ds. 

The quantities ,„(«) and — A^- ,„(n) are always positive. 

Proof The proof relies on the careful examination of the surface integrals arising in the cal- 
culations. Let e and R he a. small and a large, respectively, enough positive numbers. If 
Bfrji = Bii\Bir^rn- Calculating the norm of u in J" on -Be^j? we face the existence of the fol- 
lowing surface integrals: 

/ \x\^^ ds., I \x\~^ ds.. 



and 



jQ^J^u^ds, / Jq^ J'qU^ ds. 



For e small enough, the first surface integral corresponds to A^. The third corresponds to K^+^j^{u) 
and Aj- ,„(«). Observe that Jq = — Ji and that Jq and Ji have the opposite sign at Zm + e, 
while they have the same sign at Zm — £] thus, the quantities Aj+ ,„(«) and — A^- ^^(n) are always 
positive. 

For R large enough, the key observation is that the definition of the space J implies that 
u € L^(IR^); the second and the fourth surface integrals tend to zero. Finally, we conclude that 
(f03]) is true. ■ 

Next, we investigate some further properties of the functions that belong in J and J . Observe 
that if t> G JT" and u = T^n{v), from (j6.12p we have that \\v\\j = \\u\\j > \\u\\L^^RNy Let 
De,m = -B^^+i-E(0)/-B2^+e(0), for some fixed m. Then, 



f 

(6.14) lD^mbA>-c* — prdx > -c||ti 

Jd,.^ \x\' 



2 

2 



i.e. is bounded from below, for every m. Assume, for the moment, that u = 0, in B^{0). From 
()6.14p and (j6.13p we conclude that n G J', or equivalently u G JT", if and only if 

oo 

(6.15) K+,mi'a) - K-,m{u) < OO. 

m=l 

This actually means that 

(6.16) lim u{zm + e) = lim u^z^^i — e) = 0, 

e->-0 e->-0 
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for any m, and particular 

(6.17) limu(zi-e) = 0. 

Then u € ^.^{D^^m), for any m, and Iqc^q^[u] is a well defined positive number. 
The behavior at the origin, as it is expected, is a local effect and is described in Section 2. To 
see this we may use the standard argument using cutoff functions; /_b^(o) — is at least bounded 
from below. We finally have 

Lemma 6.3 For every J and every e, m, the quantity Ib^^[u] is well defined and positive. 
Then, u ^ J if and only if 116. 15\) holds. Moreover, u vanishes at Zm with rate 0(||a;| — Zm\'^), 
a > 1/2. In particular, the restriction of u on -621(0) is a function that belongs to H, i.e. 

J{B^^{0)) is equivalent to H{B^^{0)). 

Its behavior at the origin was studied in Section 2. Finally, u G -L^(M^), so it tends to zero at 
infinity. As it follows, 

(6.18) if uej then u€H\B';^), 
t.e. J{B-JCH\B';J. 

Compact imbeddings From the above discussion and the compact results of (|3.2p . we have 
the following compact imbeddings 

• Jr{Bz^) ^ Lr{B^_^), 2 < q < -^T^, where Jr and Lr are the radial subspaces of Jr and Lr, 
respectively. 

. J{B,,{0))^LP{n),l<p<^, 

• J ^ Lli^^), \<q< where g is a weight function belonging in 

Conclusion The space ^7 is a well defined space, through J ^ as the isometric space of J under 
the transformation T^n . In this space there exist "bad functions", concerning the behavior at 
the origin, as was discussed in Section [21 Concerning the behavior at infinity, each function 
belongs in L^(M^) and in H^{B'^_^ ). In addition, if u G must vanish at the zeroes of the Bessel 
function, as Lemma 16.31 states. 

As we can see, the transformation (16. Sh removes the singularity at infinity. However, the cost 
we have to pay, is the existence of an infinite series of singularities at Zm- Contrary to the case 
of the origin, there exists a singularity energy at these points, which is possible to increase the 
/[n]-" energy" . 

The study of the Cauchy problem (jl.ip on is not in the purpose of this work; actually this 
was done in [14]. However, we make some comments. This problem is well defined on J'. An 
important observation is that when u{x,t) is a solution of equation (jl.ip . then v satisfies the 
following associated equation 

(6.19) vt = J-2 V • J2 Vv) - V, 

with clear equivalence for x ^ 0, z^- The associated space is J', with base space X^n = 

L2((i/X,M^), dfl = J2, with 
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However, as it seems, the above problem may be not appropriate in order to study (jl.ip . First, 
the positive solutions of (j6.19p correspond, through ()6.3p . to solutions of (jl.ip which change sign. 
In addition, these solutions must vanish at z^, thus, the dynamic of (|l.ip described by (j6.19p may 
be nothing else but the dynamics of (jl.ip defined on each bounded domain Bz^^i{0)\Bz^{0), up 
to a scaling. 



7 Appendix 

7.1 Alternative proof of Proposition 12. 11 

We will give the proof of Proposition 12.11 using a different approach from this of Subsection [27 



Alternative Proof of Lemma \2.4\ Let v he a. function such that v\qq^ = and holds. We will 

find a C^(r2)-function such that 

(7.1) / [V</>- Vvpdx < e, 

Jo. 

for any e > 0. Note that since holds, we get that V e Hl{VL), so V L'^{VL). Thus V may 

be written as the sum of its radial part Vr and the nonradial part Vnr ■ If we denote the same for 
(b we have that 



Jn J \Jn J 

(7.2) |xr(^-2) \Vcl)nr " VVnr\^ dx 

The first integral in (|7.2p is equal (up to some constant) to 



[V0- Vvpdx, 



for some R > 0, large enough. Since the C^{Br) functions are dense in Hq{Bji C M^) we have 
the existence of a Cq°{Br) function such that 

1/2 



for any e > 0. It is clear that this function is a mollifier of Vr in Br C and may be chosen to 
be radial belonging in C^(il). For the second integral in (17. 2p we have that (pnr — Vnr has zero 
radial part. Then, \x\~^^~'^^^'^{(f>nr — ^Vnr) belongs in H^{Q). More precisely, 

\x\-'^''-'^\V<Pnr-'^Vnr\^dx = I" (^"^^^ (^„, - ^;,,)] 

< rlllTr(^"2)/2^ , _ N||2 

— 111-^1 yvrir n-i' ) \\ Hq (Q) '> 

where / is the Hardy functional. However, the C^(O\{0})- functions are dense in Hq{Q,), for 
instance see (j2.15p . Thus there exists a function cpnr, belonging in Co°(r2\{0}), such that 

1/2 
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for any e > 0. Is clear from the construction of Lemma 12.11 that (j)nr has zero mean value. Finally, 
setting (f) = + (j)nr, we conclude that ()7.ip holds for any e > and the proof is completed. ■ 

We define the space as the completion of the set 

{AT— 2 1 
w; w=\x\—(j), (j)£C^{n)[ 



under the norm 



From the fact that 



, JV-2 



V x\~ 



2 



dx . 



for any € C^(il), we obtain that this space is isometric with H. 

Next we give the relation between these two spaces: H and "H. Note that for every function in 
the set B, we have that ||w||^ = Thus, from Lemma 12.41 w belongs in T-L and since the 

two norms coincide we have that 

Moreover, we have the following 

Lemma 7.1 The space H is not a strict subspace ofH. 

Proof Assume the opposite. Then there exists function v G H, v ^ 0, such that 

V ±'H, 

which means that 

< v,w >-fi= 0, for any w Gl-L. 

The latter is equivalent to 

j \x\-^^-'^^Vv -V (^\x\^ (1)^ dx = 0, for any (/> G C^(f]). 

Integrating by parts and noting that the boundary terms are zero, we have that 

div (|xr(^-2) Vv) \x\^ (t)dx = Q, for any (/> G C^{^), 
or 

(7.3) j \x\-^ (kl^^^^^ div {\x\-^^~'^^ V?;)) (^dx = 0, for any (f) £ C^{n). 
However, we know that the operator L with domain 

D, = {fen: V • (lxr(^-2) g L2(|x|2-^dx, n)] . 

is self-adjoint and bijective with compact inverse. Thus 

(\x\^''~'Uiv(\x\-^''-'^S/v)) eL\n) 
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and from (j7.3p we conclude that 

(|x|(^-2)dii;(|x|-(^-2)V7;)) =0, 

almost everywhere in fi. The fact that L is bijective implies that v = Q which is a contradiction 
and the proof is completed. ■ 

As a conclusion we have that % coincides with i.e., % = %. This in terms of u, i.e. in terms 
of the spaces % and means that these two spaces coincide. 
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